Kinks in the electronic specific heat of strongly correlated systems 
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We find that the heat capacity of a strongly correlated metal presents striking changes with respect 
to Landau Fermi liquid theory. In contrast with normal metals, where the electronic specific heat is 
linear at low temperature (with a term as a leading correction) , a dynamical mean-field study of 
the correlated Hubbard model reveals a clear kink in the temperature dependence, marking a rapid 
change from a low-temperature linear behavior and a second linear regime with a reduced slope. 
Experiments on LiV204 support our findings, implying that correlated materials are more resistive 
to cooling at low T than expected from the intermediate temperature behavior. 
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If we trace from low to high temperatures the specific 
heat capacity cy — dE/dT of a solid, it provides for a 
rich variety of information. For a metal it increases Un- 
early, cv — JqT, with the prefactor 70 proportional to 
the density of the electronic states, i.e., 70 ~ N{Ef). 
This result is also valid for correlated systems that main- 
tain a normal metallic behavior. In this case we can rely 
on Landau's normal Fermi liquid (FL) theoryjlj, which 
describes the low-energy excitations of correlated (inter- 
acting) electrons as "quasiparticles" (QP) which are adi- 
abatically connected to the non-interacting electrons. As 
a result, only a QP renormalization factor ZpL needs to 
be included in comparison to non-interacting electrons so 
that Cv = JflT with — Jo/Zfl- This description is 
so universally applicable that special attention is paid to 
any deviation occurring in the vicinity of special points 
(e.g., Quantum Critical PointsQ, where the specific heat 
shows a logarithmic T-dependence). 

Turningback to the normal case, the common under- 
standing [3| is that the next electronic contribution to 
the specific heat is cubic, ^ T^. This is of the same order 
as the contribution from the lattice degrees of freedom, 
where the prefactor is given by the stiffness of the lattice 
and the mass of its ions. This makes the "lattice" pref- 
actor much larger than the electron contribution, so that 
the cubic phonon contribution is usually dominant Q . At 
higher temperature, finally, the specific heat saturates 
with a value proportional to the number of degrees of 
freedom in the system (law of Dulong and Petit). 

In this paper we show that the above described com- 
mon understanding of the low-temperature specific heat 
of a metal needs to be markedly corrected, if the move- 
ment of the electrons is strongly correlated because of 
their mutual Coulomb interaction. Our finding is based 
on numerical solution of the Hubbard model using Dy- 
namical Mean-Field Theory (DMFT), combined with a 
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FIG. 1: (Color Online) Kinks in the electronic specific heat 
(right) for U jW = 0.8 (upper panel) and 1.0 (lower panel). 
The left panels show the DMFT(ED) results for the total 
energy iJtot from which the specific heat has been obtained as 
a numerical derivative after a spline interpolation. Also shown 
are two parabolic fits (see text) valid for y <; y* (j-ed solid 
line) and T > T* (violet dotted line) and the DMFT(QMC) 
results of Ref. fTH (green dots). 



field theory formula for the specific heat calculation given 
by Abrikosov et al. Q and recent results for the energy- 
momentum dispersion relationj^ @]. 

As mentioned above, starting point of our consider- 
ation is the half-filled single band Hubbard model, the 
minimal model which describes strongly correlated elec- 
trons on a lattice. This model is solved numerically using 
DMFT 0,[i| for a semicircular DOS with bandwidth W, 
and exact diagonalization (ED) [9] as impurity solver with 
7 energy levels in the bath. 

Fig. [1] shows the total energy Etot as a function of T 
and the specific heat cy obtained through numerical dif- 
ferentiation for a ratio Coulomb interaction (U) to band- 



width (W) of U/W — 0.8 (top panels) in the tempera- 
ture range where cv{T) is monotonically increasing 
Thanks to the extremely dense temperature mesh, our re- 
sults clearly show a rapid but continuous change of slope 
(kink) of cv at T* ^ O.OISW^, a feature entirely unex- 
pected for a normal FL [3] . This kink becomes more and 
more pronounced when electronic correlations are fur- 
ther enhanced by increasing the Coulomb interaction (we 
show in the bottom panels of Fig.[T]the case of U/W — 1), 
i.e, when moving towards the metal-to-insulator (phase) 
transition. At the same time, the value of T*, where the 
kink appears, is reduced, displaying a clear relation with 
the increasing correlations which reduce Zfl- 

A proper fit to the numerical data hence needs to 
consist of two slopes (renormalization factors) jfl and 
72 instead of a single one: cy = jflT for x < T* 
and Cv = i? -|- 72T for y > with a rather sharp 
crossover in between. This has been achieved through 
fitting EtotiT) = [EtotiO) + 7flTV2] /(T - T*) + 
{Etot{T*) + BT + -t2T^/2)[l-f{T~T*)\, using a Fermi- 
function- like change f{x) = 1/(1-1- e^^) for the crossover 
at T* . Note that this fit (solid red line in left panels of Fig. 
[T]) is valid in the temperature range where cviT) is mono- 
tonically increasing, i.e., approximatively < T < 2T*. 
To assess the reliability of our impurity solver, we com- 
pared our results with precise DMFT(QMC) data [ll| 



(green dots in Fig.[Tl second row, first panel). The com- 
parison of the total energy shows an excellent agreement 
with our DMFT(ED) calculations. Notice that a direct 
observation of the kinks in DMFT(QMC) may require a 
much finer grid in the temperature regime considered. In 
our opinion, however, a first hint for a kink is already pro- 
vided by the Taylor expansion of cy (T) in Ref. .11. : The 
coefficients of the higher order terms become huge, an 
indication that the Taylor expansion is not appropriate 
as in the presence of a kink. 

To support our finite-T numerical findings, we carry 
out an analytical theory for the surprising appearance of 
kinks in correlated systems, which relies on the knowl- 
edge of the T — Green function. The analytical ap- 
proach is based on a formula by Abrikosov, Gor'kov and 
Dzyaloshinski (AGD) for the entropy of a fermionic sys- 
tem at low temperatures |3]. The AGD formula allows us 
to compute the entropy using the low-frequency behav- 
ior of the self energy S(ti;) at zero temperature, there- 
fore it connects the dynamical information (frequency 
dependence) to the thermal response (temperature de- 
pendence) . More precisely, it relates the low-temperature 
behavior of the entropy (and consequently of the specific 
heat) to the poles of the T = retarded Green function, 
which in turn follows from the the self-energy on the real 
axis. A correlated system is expected to show a kink for 
a frequency uj* << W in the T — self-energy under 
generic conditions 0]. Also experimentally, kinks in the 
angular-resolved photoemission spectrum have been ob- 
served for several strongly correlated materials such as 



cuprates [I4I, vanadates [l3| and ruthenates 14|.ll5|. 

For a normal metal, the AGD formula reproduces the 
standard Fermi-liquid result 7 = 7o/'Z^fl, with ZpL = 
[1 — In this paper we show that AGD for- 

mula also works beyond this linear Fermi-liquid regime, 
and it actually describes a kink in the specific heat at a 
temperature T*, if the proper "kinky" T = self energy 
for a correlated electron system is used. 

Let us now prove this result. The specific heat can be 
expressed via the entropy as cy (T) — T^. For a metalhc 
system at low T, the entropy is computed according to 
AGD 3 as 



S{T) 



deN{e) 



[logG],i(e,c.)-logG^i(e,c.)], 



(1) 



where N{e) is the non-interacting DOS (We use a semi- 
circular DOS N{e) = ^V-D^ - £2) with bandwidth 
W — 2D, f{uj) = Quj/T^i is the Fermi-Dirac distribu- 
tion function and G/j/^(e) the retarded/advanced T = 
Green functions respectively. Eq. ([T]) has been obtained 
in Ref. [sl by a low-temperature expansion of the Self- 
Energy. Introducing the auxiliary dimcnsionlcss variable 
y = Lo/T (ks = 1), and performing a straightforward 
derivative w.r.t. T, the specific heat is eventually com- 
puted as 



cv{T) - = / deN{e) / dyy 



:[GA{e,yT)—G^\e,yT) 



GR{e,yT)-^G'^'ie,yT)]. 



(2) 



In the case of a FL, where just one renormalization factor 
ZpL is present for the low-frequency behavior of the self- 
energy, the standard FL formula (cy(T) = jq/Zf^T) is 
easily recovered. 

The same Eq. ([2]), however, yields completely different 
results for strongly correlated metals: When the interac- 
tions are strong enough, the spectral function displays 
a typical "three-feature" structure (the QP peak, and 
the two Hubbard subbands), which survives to moderate 
doping. In this situation two distinct renormalization fac- 
tors can be identified in the low-frequency regime with 
a kink in the real part of S(w) in between [6|. Specifi- 
cally, Ref. shows that while the lowest frequencies fol- 
low the FL behavior ReE(aj) = (1 — 1/Zfl)^ there is a 
rapid (but continuous) change of slope (kink) in Rel](a;) 
at frequency w* < (e.g., lu* ~ (^/2 - 1)ZflW/2 in 
the case of the semicircular DOS). For larger frequency, 
Rel](a;) = ~b + {1 — 1/Zcp)uj with a reduced slope 
Zcp > ^FL (typically by about a factor 2), see inset of 
Fig. [Hand Ref. [6]. The constant b = (1/Zfl -1/Zcp)uj* 
ensures the continuity of T,{uj). 



As a consequence of this self energy kink, the Green 
functions and their temperature derivatives appearing 
in Eq. (0) have to be written separately for the two 
regimes, namely Gr = {yT/Z-pi^ — e + iO^)^^ and 
-^G^{e,yT) — y/ZpL for w < w*, while for frequencies 
larger than uj* one has Or = [yT/Zcv - e + 5 + iO+)~^ 
and -^G~^{e,yT) = y/Zcp with Zcp > ^fl- The pa- 
rameters have been extracted from fitting S(cci) of Ref. 
(Til (the use of the numerical renormalization group as an 
impurity solver allowing for very accurate low-frequency 
results). In the inset of Fig. [5] we show and the fit 

(blue line) . We recall in passing that it is Zqp , which con- 
trols the width of the "quasiparticle" peak in the inter- 
acting DOS; while Zfl characterizes only the asymptotic 
properties in the limit oj — s- (or T — s- 0) IS] . 

The evaluation of Eq. ([2|) has been performed by split- 
ting explicitly the integral over y in the two regions: 




This equation is the final result of our analytical calcu- 
lation. It allows us to compute, through a simple inte- 
gral, the specific heat from the non-interacting density- 
of-states N{E), the two renormalization factors Zfl and 
Zcp and the kink frequency uj*. Using the parameters 
extracted from Ref. [la| we obtain the solid line shown 
in Fig. H 

It is easy to verify that the standard Fermi-liquid be- 
havior is recovered from Eq.[3]in the limit of large lu* (i.e., 
when only one low-frequency scale is present) . In the op- 
posite limit — !■ a standard Fermi-liquid behavior is 
also recovered, though with a different renormalization 
factor 7 — jq/Zcp. More interesting is the intermediate 
situation, which we are considering here, when lu* lies in 
the low-frequency range. In this case, the specific heat 
behavior shows a kink at a temperature T* cx u* (with 
a proportionality factor of about 1/5 for the case of the 
semicircular DOS): As one can see in Fig. [2]the standard 
FL behavior cy(r) = is recovered only for T <T* , 

while at T = T* a sharp change of slope is observed. 
For T > T* , the specific heat is still essentially linear, 
but with a completely different slope (determined by the 
value of Zcp and the coefficient b). Eventually, the AGD 
formula loses its validity at higher temperatures, where 
the maximum of cy is reached (see again [lo|V 

We emphasize that the AGD formula reproduces the 
finite low-T DMFT(ED) solution not only qualitatively 
but also at a quantitative level. This allows to precisely 
relate the value of T* with lu* and, hence, with the char- 
acteristic parameters of the system (e.g., the estimate 
T* - l/10(\/2 - 1)ZflW works well for the case of a 
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FIG. 2; (Color Online) Analytical theory describing kinks in 
the specific heat (red solid line) on the basis of the AGD for- 
mula (see text). The blue dashed line in the inset was fitted 
to the (red) numerical renormalization group data of Ref. [l6l : 
note that the deviation at larger frequencies uj does not sig- 
nificantly affect the specific heat in the plotted temperature 
range. The agreement of the analytical calculation with our 
numerical results (black crosses) is excellent. 

semicircular DOS). The agreement with the AGD for- 
mula is particularly remarkable if we notice that, strictly 
speaking, the AGD formula is only applicable to the 
linear-T regime, since it does not include all additional 
terms leading to the aforementioned T'^ contribution. 
However, this term is small at low temperatures. Hence, 
if the correlation is strong enough, it can push the kink 
in the very small T regime, where the AGD formula is 
expected to work. This explains why our analytical cal- 
culation is able to reproduce our numerical results to a 
very good accuracy in Fig. [2l Let us emphasize that it 
was not at all clear a priori whether an AGD-likc calcula- 
tion was possible beyond the regime of Landau's QP, i.e., 
after the kink in the energy-momentum dispersion which 
indicates the basic excitations are no longer Landau QP. 

The theoretical evidence of a low-temperature kink in 
the electronic specific heat of strongly correlated systems 
poses the question of its experimental observation, which 
was -so far- still lacking. The main problem is obviously 
the phonon contribution cy which -because of its 

large prefactor- usually overshadows the much smaller 
electronic contribution to the specific heat, already at 
temperatures of few ten Kelvin. This restricts the choice 
to materials which show the kink at a very low T* . That 
means in turn compounds with a strong renormaliza- 
tion (Zfl ^ l)j i-e., heavy Fermion systems. Given our 
starting point, the Hubbard model, the ideal material is 
LiV204, the first d-electron system where heavy Fermion 
behavior was found Indeed, recent LDA+DMFT cal- 
culations [18], which take into account the realistic three- 
d-band structure of LiV2 04, have demonstrated that an 
effective description in terms of the single-band Hubbard 
model (very close to half-filling) is particularly appropri- 
ate for this compound. 

In Fig. [21 we show that our theory nicely describes the 
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FIG. 3: (Color Online) Kink in the low temperature specific 
heat of LiV2 04 (open blue circles) visible at T* ~ 5 — 6K, 
and well reproduced by our analytical theory (red solid line) 

experimental results for LiV204. We compare the data 
of Ref. [3 (displayed in a magnified low-T range with 
respect to the original publication) with our analytical 
formula, fitting the free parameters to the experimental 
data. Indeed, a kink is clearly visible, as at the curve 
rapidly changes its slope at a temperature T* of 5 — 6K. 
The three fitting parameters (Zfl = 0.054; Zcp = 0.092, 
Lu* = 0.0035Vl^ with W = 600 meV for LiV204) assume 
very reasonable values. This clearly confirms the strong- 
correlation origin of the kink in the specific heat of this 
material. 

We notice there are also kinks in the specific heat 
of /-electron heavy Fermions such as YbRh2Si2 20] or 
YbCu5_a;Ala;[2H. However these materials are close to 
a quantum critical point, at which additional physical 
processes become important. In some systems also long 
range magnetic order leads to additional structures in the 
specific heat. Hence, at present, it is less clear in how far 
these kinks are connected to our theory. Another material 
with strongly correlated Fermions showing similar kinks 
in the specific heat is "^He (Ref. 12^ [23I) for which however 
the application of a lattice model such as the Hubbard 
model represents certainly quite a crude approximation. 

In conclusion, we have demonstrated numerically, ana- 
lytically and experimentally that the textbook knowledge 
of the electronic specific heat at low temperatures needs 
to be modified for strongly correlated electrons. In the 
proximity of the Mott transition the leading correction 
to the linear Fermi-liquid temperature behavior is a quite 
rapid change of slope, i.e., a kink, which takes place well 
before (at smaller T) the standard T"^ behavior becomes 
relevant. Let us emphasize the reported kink is a generic 
feature of strongly correlated electron systems, in very 
contrast to existing theories for kinks stemming from the 
coupling to (potentially present) bosonic degrees of free- 
dom . Since the slope of the specific heat is reduced after 
the kink, the behavior of cv/T is just opposite to what 
one would expect from the standard theory, i.e., cv/T 



is decreasing with increasing temperature instead of the 
expected increase due to the cubic term. Hence, if one 
extrapolates from the behavior at intermediate tempera- 
tures (i.e., after the kink) without taking into account the 
kink, a much lower specific heat at low temperatures is 
obtained with respect to the actual result. In other words, 
a material with strongly correlated electrons can be un- 
expectedly resistant against cooling at low temperatures. 
Moreover, depending on the temperature range consid- 
ered in the experiments, only one of the two regimes of 
linear behavior of cv{T) may be accessible. This can eas- 
ily lead to remarkable inconsistencies in the analysis of 
the experimental data for strongly correlated materials. 
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